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We consider the hydrodynamics of strongly interacting quark gluon plasma in 
finite temperature and density using the holographic duality of charged black hole in 
anti DeSitter space. We calculate the transport coefficients at arbitrary energy scale 
by considering the holographic screen at finite radial position. We first calculate 
the flow of sound velocity in this method and check the consistence with previous 
result. Then we calculate diffusion constant of charge and find that Einstein relation 
between susceptibility, conductivity and diffusion constants will hold at arbitrary 
slice. 

I. INTRODUCTION 

In usual application of AdS/CFT [JJ, physical quantities are calculated at the asymptotic 
AdS boundary. However, a single asymptotical AdS gravity system is expected to describe 
a theory both at the UV regime as well as the IR region, since the radial direction is 
identified as the energy scale [2HI] of the gauge theory. While the asymptotic boundary 
describes the UV fixed point, the horizon, as IR boundary, should be able to characterize 
low energy behavior of the theory. In fact, some of the transport coefficients were calculated 
using membrane paradigm at the horizon. Therefore it would be very interesting to find the 
interpolating functions connecting the UV and the IR behavior. Such interpolating functions 
can be computed by putting the holographic screen at a slice r = r c . The paradigm to 
compute the slice dependent quantities has been recently proposed in [51 E] as holographic 
Wilsonian renormalization group. Physical quantities at the slice have been defined in a 
slight different way of Wilson fluid/gravity duality in [7]. It has been shown that these two 
different methods, together with the earlier sliding membrane paradigm [8] are consistent [9j 
ITU] . For further recent progress in this direction, see [ITl - fTT] . 

When we ask the density dependence, we need to introduce the local U(l) charge, which 
in turn gives an extra complexity in the hydrodynamic analysis due to the mode mixing: 
in the RN-AdS black hole, vector modes of gravitational fluctuations mix with transverse 
Maxwell modes [21] [23]. The scalar modes of metric perturbation mix with longitudinal 
modes of Maxwell potential and the analysis is much more involved [22]. In a companion 
paper [10], we focused on the former, i.e, the holographic RG flow of tensor and vector part 
of metric perturbations and transverse Maxwell perturbation. We proposed a systematical 
method to compute the hydrodynamical poles and retarded Green function at the slice 
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In this paper, we focus on the running of the sound modes. Notice that the running of 
sound modes has already been studied in [7j and in more recent paper [13]. What is new 
in this paper is the density effect on the sound modes coming from the mixing of gravity 
and Maxwell part in the presence of the charge. We will obtain analytical results for sound 
velocity and charge-diffusion constant and will see that they have nontrivial charge depen- 
dence in RG flows. We also find that Einstein relation between susceptibility, conductivity 
and diffusion constants will hold at arbitrary slice. 

This paper is organized as follows. In section [TTJ we analyze the sound part perturbations 
for the RN-AdSs and focus on the dynamical structure of Einstein and Maxwell equations. 
In section III, we write down the master equations which are equivalent to equations of 
motion in section [TT] and solve them perturbatively in hydrodynamical limit. We impose 
both the in going horizon condition and Dirichlet boundary conditions at the slice surface 
r = r c . In the following section IV we derive the on shell action for the sound fluctuations 
and write down forms of Green functions and pick up the pole position. We found a charge 
diffusion pole and a sound pole at r = r c and we discuss the RG flow of them. We found the 
universal value (charge independent) located at the horizon for both momentum diffusion 
and charge diffusion constant while located at the UV boundary for the speed of sound. We 
checked some transport coefficients and Einstein relations at the slice in the last section. 



II. SOUND MODES IN REISSNNER-NORDSTROM-ADS 5 

We consider the 5 dimensional Einstein-Maxwell theory. The action is given by 

S = ^ / <fx^g (R-2A)-±J a*xy^g F 2 , (1) 

where five dimensional Newton constant is given by G§ = k 2 /87t. The metric of the RN-AdS 
black hole is given by 

ds» = r ^(-f(r)dt* + dx 2 ) + ^, (2) 

where 

/M = l-^ + ^, A = -«, (3) 

and the background gauge potential is 



r 

where 

e 2 2k 2 ' 

The chemical potential /i is fixed by the regularity condition at the horizon as 



(5) 



H = -J. (6) 

' 4- 
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It is convenient to introduce a coordinate u = r^/r 2 , where r + is the outer horizon. Then 
the metric is expressed as 

I 2 l 2 dv 2 
ds 2 = -L- {-f(u)dt 2 + dx 2 ) + (7) 
Ab 2 u v v ' 1 4u 2 f(u) 

f(u) = (1 -u)(l +u-au 2 ), (8) 

where 

l 2 q 2 I 2 



(9) 



Consider the gravitational fluctuation for the RN-AdS solution, the Gibbons-Hawking term 
on the boundary is necessary to well define a variation principle 

S GH = ^[d A x^lK (10) 



K 2 j 

where 7 M „ and K^ v are induced metric and extrinsic curvature on the boundary, respectively 
Furthermore, in order to regularize the on shell action, we introduce the counter term [20J 



S« = ±J*xV=lj. (11) 

From now on, we study fluctuations h mn , B m for the background metric g mn ^ and the 
gauge field A 

9mn 9mn h mni A m -A m -\- , (1*^) 

[I 

and consider the linearized theory of these fluctuations. Take the gauge condition h um = 
and B u = and consider the Fourier transforms of these fluctuation fields, 

M*) = / ^i e ~ MkZh ^) > ( 13 ) 

where the momenta k is along the z-direction. The fluctuations of the metric are categorized 
into the tensor mode, shear mode and sound mode, and the gauge field is categorized into 
the vector mode and scalar mode. At the linearized order, there are mixture between the 
shear part of metric and transverse gauge field, also the sound mode of metric and the 
longitudinal gauge field. Here, we focus on the sound mode of the metric and the vector 
mode of the gauge field. Then, the relevant components are 

hi, hl = hl h*„ hi B u B z (15) 

where we have defined h^ = g^ x h\ v for later convenience. 
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Relevant equations for the sound modes in the Einstein equation are (t,t), (t,u), (x,x), 
(z,z), (u,u), (t,z) and (u,z) components listed as follows: 



* U f V 2 2 / U f 



2b 2 ( 1 \ 7/ 7/ 

+ f w 2 /i* + -A z z + wfc^J + 2a-/i* + 4ajB' t , (16a) 

o = K" + ^ "^V - ^ (tf + tf) + £ K - *V) 

- - 2a^, (16b) 

o = ftf + ^to^r - - (\ht + hA + ^ + 2 ^ 2 - k^fhi - 2efK) 

U 2 J M \Z J Uf 

11 11 

- a-h\ - 2a^, (16c) 

= fcf - -/if + - 3auB' (16d) 

it it/ 

= h? + 2hf + K" + f j (^ht + h%' + \h'j} + 2a^ + Aa- f B[, (16e) 



= w 



2 ^' + *r - j (k + ^: 



+ fc ( h? - f -hl ) , (16f) 



= kh\ + 2kh x J - 7 /if + -y^ + 3a- (kB t + wE 2 ) . (16g) 

The Maxwell equations which is relevant to the vector mode are t, z, and u components: 

= B'l - ^ (k 2 B t + kcoB z ) + l - (hi' - 2h* x ' - h?) , (17a) 

= B';+ f jB' z + ^- {u 2 B z + ukB t ) - -hf, (17b) 

= uB' t + kfB' z + ^(hl- 2h* - K) - kh z t . (17c) 

There are 10 equations for 6 fields, and only 6 of these equations are independent. These 
equations are classified into the dynamical equations which are second order differential 
equations and first order constraints. For the Einstein equations, we have 4 dynamical 
equations and 3 constraint equations. Independent equations are 1 dynamical equation and 
3 constraint equations. For the Maxwell equations, 1 dynamical equation and 1 constraint 
equation are independent. These differential equations give 8 integration constants since 
the solutions of n-th differential equations have n integration constants. We impose the 
incoming boundary condition at the horizon for each dynamical equation. Then, there 
remain 6 integration constants which are determined by choosing the values of the 6 fields 
on the UV boundary. 
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III. MASTER EQUATIONS AND BOUNDARY CONDITIONS 



A. Master equations 

Consider the two independent dynamical equations. By using an appropriate redefinition 
of the fields, we obtain the dynamical equations which always satisfy the constraint equa- 
tions. Such fields and equations are called master fields and master equations, respectively. 
The master fields are derived in [2T], and can be expressed in our notation as 
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-[(Ab 2 k 2 -3f)K + 2f {K' + K')] 



Av?l 2 (Ab 2 k 2 -3f) 
for metric perturbations and 

A = 2a (-hi + 3h x x - 2B' t ) 



(18) 



(19) 



for gauge part. Two master equations for the master fields ( 18 ) and ( 19 ) can be diagonalized 
by the following new master fields: 



where 



u 1/4 

$ ± = a± $ + -—-A, 



a± — (1 + a)(l ± a) — 3au c 



Aab 2 k 2 
(1 + a) 



(20) 

(21) 
(22) 



The final master equations become 



(u^ 2 fY 



(23) 



where potentials V± are given by 



lQu 2 f 2 (4b 2 k 2 -3f) 2 



x | - 4uf(4b 2 k 2 - 3f) (\Qb 2 k 2 {b 2 k 2 - C±u + 3au 2 



- Af(2b 2 k 2 + 3C±u - 9au 2 ) - 3{f) 2 
+ f 2 16 (b 4 k A + l2C±b 2 k 2 u - 108ab 2 k 2 u 2 + 54C±au 3 - lQ2a 2 u 4 

b 2 u 2 



24/' b k 1 - 6C±u - 18au 2 + 9(/ 



p/\2 



+ 



Uf 2 



(24) 
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Let us solve these master equations. We first impose the incoming boundary condition 
at the horizon. Near the horizon u = 1, the solutions of the master equation (23) behave as 

$± ~ (i - u y, (25) 

where there are two independent solutions with v = ±ibu/(2 — a). Each of them corresponds 
to the incoming and outgoing modes at the horizon. We impose the incoming boundary 
condition by taking only v = —ibu)/(2 — a). Then, the solution can be factorized as 

$ ± = (1 - u )- ib " /{2 - a) F±. (26) 

Since we have factorized the most singular part near the horizon, F± should be regular at 
the horizon. In the hydrodynamic region, u and k are small and we can solve the master 
equations order by order by expanding the master fields with respect to u and k. It was 
solved in [22], and the solutions take the form of 



$ ± = C±H±{u){\ - u )- lhui/{2 - a) (1 + buF± w (u) + b 2 u 2 F± 20 {u) + b 2 k 2 F± 02 {u) + •••). (27) 
where H± given by 

H + = u~y\ H = , (28) 

(1 + a) — ^au 

are factorized such that the leading terms in F± become u independent constant . The deter- 
mination of integration constants C± is important for computing the final Green functions 
and we will see it soon. Details of F±ij are shown in Appendix A. 



B. Boundary Conditions at Slice u = u c 

In order to determine the integration constant C± we need boundary conditions of the 
original fields and B^. Consider the boundary at slice u = u c and impose the boundary 
values for the original fields as 

K(u c ) = h^°\ B^u c ) = B^. (29) 

Taking use of the Einstein equations and Maxwell equations at the slice, we obtain the 
following relations between the integration constants and boundary conditions: 

£ + C + + £_C_ = £o, (30) 
M+CV + M_C_ = M , (31) 

where So and A4o are linear combinations of h^ ^ and (B^)^ . Namely they take the form 
of 

£o = £oi<pf\ Mo = M Kpf\ (32) 
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where <pi stands for the original fields, h% or B^, and £qi,Moi are given functions of uj, k, 
a, b and u c . £± and Ai± are also known functions of uj, k, a, b, u c and generally they take 
the form of 



£± = £± F±(u c ) + £±iF' ± {u c ), 
M± = M± F±{u c ) + M±iF' ± {u c ), 



(33) 
(34) 



where £± , £±±, A4± , Ai±i are all known coefficients and F± are defined in terms of the 
solutions of the master fields as 

F±(u) = (1 - u)-^/^ (1 + buF ±w (u) + b 2 u 2 F ±20 (u) + b 2 k 2 F± 02 (u) + •••). (35) 

Remember F±(u) are given functions since we have solved the master equations. Once the 
above ingredients are given, C± can be finally expressed in terms of the boundary conditions 

as 



C 4 



£ M- - M £- 
£ + M--M+£- 



C- = 



M £+ - £ M+ 
£+M--M+£- 



(36) 



IV. POLE OF GREEN FUNCTIONS AT THE SLICE u = u c 



Now we consider Green functions. By using the equations of motions, the on-shell action 
for fluctuations is reduced to surface terms as 



_J! f d * k J -h z h zl + Lh x h x ' + f-tfh*' + —h l h zl 

32/^4 J (27r) 4 I u + u + u h ^ + 2u thz 



+ f -h x X + f -h x x h z J + ^-h'X + f -h z z h x x 
u u 2u u 



(f - Vl) {hi? - A (6/ - uf - 6 y/f) hX 



u 2 f 

1 

4m 2 



l - 2 (6f - uf - 6y/f) KK - 1 (3/ - u/' - 3 v7) ^ 

) 



+ Za(B t B' t - fB z B' z + l -B t h\ + B z h\ - B t h x x 



BtK 



u=u c 



It can be rewritten as 



/dk 4 ( \ 
j—^ ^ 7 (u c )/C 7J 0j(w c ) + <Mit c )£jJ0j(«c)J • 



(37) 



(38) 



By using the constraint equations and the definitions of the master fields, we obtain the 
relation between the original fields and its first order derivatives 0' 7 . Since we have 4 
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constraint equations and 2 master fields, we can diagonalize these equations such that the 6 
first derivatives 0' 7 at the boundary can be expressed in terms of the boundary values 
They take the following form: 



*>i = n IJ( j>j + n I+ c + F + + n^c.F. 



Combined with (38) the on shell action is obtained as 

dk A 



S 



(2tt)< 



J > 



where Qu is given by 



ICu + £ik I U K j + U K +F + (u c ) 



SqjM- - Mqj£- 

S+M- - M+E- 

+ Hk-F- [u c ) - — — - 

£ + M- - M+£- 

Finally the retarded Green functions at the slice G R are given by [19] 

Gfj = 2Q U . 

We find the pole of the Green functions is located at 

E + M- - = 0. 



For convenience we rewrite (35) as 



F±(u) = 1 + bcoF ±10 (u) + b 2 u 2 F ±20 (u) + b 2 k 2 F ±02 (u) + 
The pole has the following structure: 

= P 40 w 4 + P 3 ou 3 + P22to 2 k 2 + +P 12 cok 2 + P 0i k 4 + ■ ■ 
where the coefficients are listed as follows 

P 4 o = 9b 2 f'(u c ) [f +10 F'_ 10 + F'_ w (u c ) 
P30 = 9bf(u c )F'_M , 



22 



-3b 2 [f(u c )(f'(u c ) + 6cm 2 ) - u c f 2 (u c )} 



F +10 F'_ 10 



F'-20{U C) 



54a 2 b 2 u 2 c 



(39) 



(40) 



(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
(47) 
(48) 



+9b 2 f'(u c )F'__ 02 (u) - 12b 2 f 2 (u c )F' +l0 (u c )F'_ 10 (u c ) + ^ + (/)/7(// 
P12 = -36 [f(u c )(f(u c ) + 6m/ 2 ) - u c f 2 (u c )] F'_ 10 (u c ) , 
P 4 = {-3b 2 [f(u c )(f'(u c ) + 6au 2 c ) ~ u c f' 2 (u c )] FL 02 (u c ) + {1 + a)f/2 

6a(l + a) + (8 + 24a - 3a 2 + 8a 3 )u c - 18a(l + a) 2 u\ + 15a 2 (l + a)u\ } . (52) 



x 



(49) 
(50) 
(51) 
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Apparently Pij coefficients here can only be fixed up to an overall constant and we leave the 



explicit form of S+Ai- —S-A4+ to the Appendix B. Due to the structure (45), the pole has 
4 roots which are referred to as 



cj — 1/r ^ , u = ±c s k + • • • , 

and can be factorized to be product of these roots: 



-iD A k 2 + 



= (u - c s k H )(uj + c s k H )(u + iD A k 2 H )(w - 1/r H ) 

= - l -^ + u 4 + c ^ujk 2 - (c 2 s -i^-)u 2 e-i^ k A + --- . 

T T V T } T 



(53) 



(54) 
(55) 



Here "• • ■ " denote the higher derivative corrections. By matching (45) and (55), we obtain 
the bare sound speed 

2 _Pi2 _ f(u c )[f'( u c ) + 6aul] - uj ,2 (u c ) 

C s p 
r 30 



3/'(« c ) 



(56) 



and the bare charge diffusion constant 
.-Po4 



-i— — = Im 

-Pi 2 



} FL 02 (u c ) 



2ong[6o(l + a) + (8 + 24a - 3a 2 + 8a 3 )w c - 18a(l + a) 2 M 2 + 15a 2 (l + a)u 
' (1 + a)f'\u c ) [f(u c )(f'(u c ) + Qaul) - u c f'\u c )} F'_ w (u c ) 



(57) 



In the orthonormal frame, the proper frequency u c and proper momentum k c are defined by 

oj . k 



V~9tt ' y/9zz 

The normalized speed of sound and the diffusion constant are given by 



_ 2 _ _9zz 2 _ f( u c)[ f'(u c ) + Qau 2 ] - u c f' 2 (u c ) 



D f 



9tt 

Qzz 
9tt 



T H D A = Im 



3/K)fK) 

4tt f( Uc )FL w (u c 



(58) 



(59) 



o(2 - a)u^[6o(l + a) + (8 + 24a - 3a 2 + 8a 3 )w c - 18a(l + a) 2 u 2 c + 15a 2 (l + a)u 



2tt(1 + a)f(u e )f'*(u c ) [f(u c )(f'(u c ) + 6au 2 c ) - uj'*(u c )] F'_ 1Q {u c 



(60) 



V. RG FLOWS OF SOUND SPEED AND DIFFUSION CONSTANTS 



We will summarize radial flows of various quantities in this section. It is shown in Left 
part of Fig [T] that the bare speed of sound flows quite differently for different charges. For 



10 



zero charge case, c 2 s {u) runs from | to | if we move from UV boundary to IR horizon. 
This is consistent with the chargeless result in [3 [T3j at the horizon and also at the UV 
boundary [IB] . As suggested in [7], the physical speed of sound may be defined as the one 
in the orthonormal frame, which is shown in right part of Fig [TJ The normalized sound 
speed always diverges at the horizon for any charge. One interesting observation is that 
both for bare one or normalized one, there is a universal UV boundary value of sound speed 
4g for different charges. Let us turn to the diffusion. Remember the momentum diffusion 
Dh is encoded in hydrodynamical pole in shear part of metric fluctuations [10J and charge 
diffusion Da is contained in sound part. As shown in Fig [2J both for normalized Dh and 
Da, the universal charge independent value 4- is located at the horizon. 




FIG. 1: Radial flow of speed of sound (? s . Left: The bare speed of sound. Right: Normalized 
speed of sound in orthonormal frame. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

LI U 



FIG. 2: Radial flow of diffusion constants. Left: The momentum diffusion constant coming from 
shear part hydrodynamical pole [10 i. Right: Charge diffusion constant coming from sound part 
hydrodynamical pole. 
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VI. TRANSPORTS AND EINSTEIN RELATION AT THE SLICE r r 



Retarded Green functions can be computed explicitly following (41 ). Generally the Green 
functions have complicated expressions. Rather than showing all the explicit expressions we 
focus on G zz , Gujj and Gtt at certain limits. By using the Kubo formula, we obtain 

a = lim —lmG zz (u, k — 0) 



u-+auj ■ 16k 2 6V 2u c f{u c )\mF'_ 



10\ u c 



= l —, . (61) 

2e 2 bu c f (u c )lm.F'_ 10 (u c ) 

This is electrical conductivity, which is the same as our previous result \±0\. We also find 
that the bulk viscosity vanishes at arbitrary slice u = u c 

C = lim - ImGiijj(u, k = 0) = . (62) 

i ,j=x,y,z 

In k — > limit, the density-density correlation function Gtt behaves as 

r lmG tt __ 13 3a6 

fc™ k 2 " WK 2 fi 2 2cou c f(u c )lmFL w (u c )' 

Together with the pole structure which we have discussed above, it implies that 

tt[UJ ' ] ~ 16&*kV 2u c f(u c )lmF'_ w (u c ) (^ + D\V){u? - c 2 k 2 ) + " ' [M) 
and then we obtain 

ImGUw, k = 0) = lim lmG tt (u, k) = J— ^4 ^-S(u) . (65) 



By using the definition of the charge susceptibility 

,_, If duj lm.Gtt{^i k = 0) 



T ./ 2vr e^/ T - 1 

we obtain 

1 / 



£>a 2e 2 bu c f(u c )\mF'_ 



10\ u c 



which satisfies the Einstein relation 



(66) 



(67) 



(68) 



D A 

As a summary, for the transport coefficients, the shear viscosity keeps a constant from the 
horizon to the UV boundary and DC electric conductivity has a nontrivial radial flow as 
explicitly shown in [10J, . Here we checked that the bulk viscosity vanishes at arbitrary slice 
and the Einstein relation holds at the slice. 
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Appendix A: Concrete expressions 

The solutions of F± are the following: 



F +W (u) = ^-^ (log (1 + u - au 2 ) - j , (Ala) 

FMu) = U4^L~ 7 1 , ), (Alb) 

+ V ; 3 \^TT4^ (l + a)uj V ; 

i-u 1 



(1 — u)(l + u — au 2 ) 

(1 - u)(l + au) log (1 + u - au 2 ) 3(1 - u)(l + au)K 2 (0) 
2(2 -a) 2 2(2 - a) 2 v / TT4^ 



x {l-u+ 



(l + a)K 2 (l)u 



(3 + (5 + 3a - 6a 2 + 2a 3 )u - 3au 2 ^)K 2 (u) 



VI + 4a 2(2 - a) VI + 4a 



(Ale) 



i^-ioH = 2 (2-a) 2 { 8(1 + a)2 l0g(M) " (2 + a)(1 + 4a) bg { l + u - au2 ) 

- 2y/l + 4a(2 + 5a) A^! (u) } , (Aid) 
3a 2 w 2(l + a)(2 + a)log(u) 



2(l + a) 2 (l + a-^a«) 12 a) 

(l + a)(2 + a) log (l + M - aM 2 ) 2 (2 + 5a + 6a 2 )A 1 ( M ) 
+ (2 - a) 2 + (2-a) 2 v / TT4^ 

i 7 -2o(«) = I du 



(Ale) 



x 



2(2 - a) 4 (l + 4a) 3 / 2 (l -u)u(l + u- au 2 ) 
8(2 - a)(l + a) 2 (l + 4a) 3/2 w(l +u- au 2 ) log(u) 

- (2 - a)(l + 4a) 3/2 (4(1 + a) 2 + (2 - 3a - 8a 2 )w + (2 + 9a + 13a 2 )w 2 

- a(2 + a)(l + 4a)w 3 ) log (l + u- aw 2 ) 
+ (1 + 4a) 2 (2 + 5a)A 2 (0)(l - u) (4(1 + a) 2 + (2 - 3a - 8a 2 )w - (2 - a)aw 2 ) 

- 2(1 + a) (2 - 2a + 41a 2 )A 2 (l) (2(1 + a) - 3aw) 2 

- (2 - a) (a(l + 4a) 2 (2 + 5a)w 3 - (2 - a)(l + 11a + 46a 2 + 18a 3 )w 2 

- (2 + 9a + 180a 2 + 224a 3 + 24a 4 ) u 

- 4(1 + a) 2 (l - 10a - 2a 2 )) K 2 (u) I. (Alf) 
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where 



2au 



K l{ u) = - Iog(l + u - an 2 ) - log ^1 - ) 
Then, F± are related to these solutions as 



F±w(u) = F±io - ; 
^±20 («) = F± 20 (u) 



log(l 



log(l - u)F±io(u) 



1 



2(2 -a) : 



log 2 (l- M ) 



F±o2(u) = F± 02 (u) 
Sqi and A^o/are defined as coefficients of the original fields in 



£o — Zoi^P , 



Mo = Moi(j>f ) . 



They can be read off from the expressions of £ and Mq: 



S = 



- 18au c (kB t {0) - ku c h x J 0) + ^ (0) ) + 3kf'(u c )(-hf ] + h x J 0) ) 



M 



+ u c {Ab 2 k 2 - 3f'(u c )) 2 2ukh z t {0) - (cu 2 + k 2 u c f'(u c ))hl i0) + cu 2 h 



kf{u c ) \-18au 2 c (kB t (0) - ku c h x J 0) + w5 2 (0) ) 



2^(0) 



(A2) 
(A3) 



(A4) 

(A5) 
(A6) 

(A7) 



(A8) 



8kWf 2 (u, 

- Ab 2 k 2 (2kB t ^ - ku c hf ] - ku c h x{0) + 2uB z W) + 3ku c f(u c )(h x J 0) - hf ] ) 



+ 



u{Ab 2 k 2 - 3f{u c )) \-2ukhf 0) + (u 2 + k 2 u c f\u c ))h x x {0) - cu 2 hf ] ] X . (A9) 



£± and Ai± contain F±(u c ) and their derivatives. They are referred to as £± , £±i, Ai±o 
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and .Mil, respecctively, and obtained as 
1 



£+i 
£- 



£-i 



M +1 
M- 



M-i 



12(1 + a)au c f 3 (u c )k 2 



-3k 2 uJ(u c )(4b 2 k 2 - 3f(u c ))(4b 2 k 2 + 2u c a_ - f'(u c )) 
+ u c (4b 2 k 2 - 3/V)) W + u c f'(u c )k 2 ) 



12k 2 f 2 (u c )(4b 2 k 2 + 3(1 + a)(l - a)u c ) 



Ab 2 k 2 - 3f{u c ) 
(1 + a)af(u c ) 



(A10) 
(All) 



x 



Q(l + a)ak 2 P(u c )f' 2 (u c ) 

x ^2k 2 f' 2 (u c ) (3u c f(u c )(3a+u c + 8b 2 k 2 ) - I8f 2 (u c )) + Ab 2 u c (2b 2 k 4 u c - 9co 2 ) 

- \2k 2 f\u c ) \2b 2 k 2 uJ(u c )(a + u c + 2b 2 k 2 ) + f(u c ) (3a + u c - 4b 2 k 2 ) - 4b 4 k 2 u 2 u c 

- 48b 2 k 4 u c f 2 (u c )f"(u c ) - 3uJ' 3 (u c ) (8b 2 k A u c + 3k 2 f(u c ) - 9a; 2 ) + 9k 2 u 2 J'\u c ) 

(A12) 

(A13) 



2u 2 c (4b 2 k 2 -3f'(u c )) 



(1 + a)af(u c )f'(u c ) 



4a(l + a)au 2 c k 2 f 2 (u c ) 



A(l + a)(l-a)f 2 (u c )k 2 

+ au 2 J(u c )k 2 (-4b 2 k 2 - Qu c a_ + 3/'(u c )) 

+ au 2 c (4b 2 k 2 - 3/K))(3a; 2 + u c /'K)A; 2 ) 



a(l + a)u c a 



(A14) 
(A15) 



2a(l + a)aA; 2 / /2 (M c )/ 2 K 



1 l '-4A; 2 a + M c / /, (M c )/ 2 (M c ) 

+ 2k 2 f'(u c ) -6au 2 c b 2 u 2 + au\(2b 2 k 2 + 3u c a+)f{u c ) + 2(a + - 3aw c )/ 2 (M c ) 

- au 2 J l2 (u c ) [4u c b 2 k 4 - 9a; 2 + 3A; 2 (/( Mc ) - u c f{u c ))] 
2a + u c 



a(l + a)af'(u c ) 
The first derivatives of original fields are expressed as 

<t>'i = nut j + n I+ c + F + + Ui-C-F-, 



(A16) 
(A17) 

(A18) 
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whose explicit forms are as follows: 

It , 2b 2 k 2 lt 2b 2 u 2 iz Ab 2 ku iz 



c/2 



f (-18au 2 + 4b 2 k 2 - 9/') - 4b 2 k 2 uf + 8b 2 u 2 + 3uf 
+ K 

3/ (Ab 2 k 2 - 3/' + 2ua.) + uf ( 3f - Ab 2 k 2 ) r< ~ 



3(a + l)a/ 2 

2u 2 (3/ {4b 2 k 2 - 3f + 2ua + ) + uf (3f - Ab 2 k 2 )) 



3(a + l)af 2 f 
Ab 2 k 2 -3f u 2 (3f'-Ab 2 k 2 ) 3uf'-Ab 2 k 2 u 

X ~ 12 / * + 12/ 2 fc 2 n * + Qf2 k 



C-F_ (A19) 



/A; 2 (18aw 2 - Ab 2 k 2 + 9/') + (46 2 A; 2 - 3/') {k 2 uf + u 2 ) _ 3au _ 3auu 
+ Upk 2 x ~ ~2j 1 " ~2fk z 



(Ab 2 k 2 - 3f) (k 2 uf + 3a; 2 ) - 3/A; 2 (46 2 fc 2 - 3f + 2a.u) ~ 
6(a + l)a/ 2 A; 2 + + 

u 2 ({Ab 2 k 2 - 3f) {k 2 uf + 3u 2 ) - 3fk 2 [Ab 2 k 2 - 3f + 2a + u)) 



3(a + l)af 2 k 2 f> 
_ 3f'-Ab 2 k 2 u 2 (Ab 2 k 2 -3f) u(Ab 2 k 2 -3f) 
6/ t+ 6f 2 k 2 z + ~ 3f 2 k 



C-F_ (A20) 



2f (9ak 2 u 2 + Ab 2 k 4 ) + (Ab 2 k 2 - 3f) (k 2 uf + u 2 ) 3au 3auu 

epk 2 K + ~r Bt + ~jk~ Bz 

(Ab 2 k 2 - 3f ) (k 2 uf + 3u 2 ) - Qa„fk 2 u ~ 
+ 3(a+l)a/ 2 A; 2 + + 

2u 2 ((4fe 2 A; 2 - 3/Q (k 2 uf + 3c 2 ) - 6a + fk 2 u) ~ 

+ 3(a + l)a/ 2 A; 2 f C ~ F - (A21) 

a;(4W-/Q a£ f 
ftt 2/ife * + 2//T 2+ / ' 

_ 46 2 A; 2 a; - 3uf ~ _ 2u 2 u (Ab 2 k 2 - 3/Q ~ 

aafk + afk + + (a + l)afkf [ ) 

B' = --h\ + -/£ + . a ~ + - 2m2q + - (A23) 

* 2*2* a(a + l)au + + a(a + l)af V ; 

K = ~ijk K + 2^Jk K + 1 K ' a(a + l)afku C+P+ ~ a(a + l)^kf f^' im 
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Appendix B: Pole Structure 



The pole is expanded as 

9bf(u c )FL w (u c )cu 3 + 9b 2 f'(u c ) 



a(l + a)P(u c ) 



F +W F'_ 10 + F'_ 20 (u c ) 



CO 



- 3b [f (u c )(f (u c ) + 6au 2 c ) - u c f' 2 (u c )] F'_ w (u c )uk 2 



-3b 2 [f(u c )(f(u c ) + 6au 2 c ) - u c f' 2 (u c )} 



F +10 FL 10 + F'_ 20 (u c ) 



10\ U c 



+ 9b 2 f(u c )F'_ 02 (u) - 12b 2 f 2 (u c )F' +10 (u c )F'_ 
54a 2 b 2 u 2 c 1 22 
+ (l + a)/'Wr" 
+ {-3b 2 [f(u c )(f'(u c ) + 6au 2 c ) - uj' 2 (u c )} FL 02 (u c ) 



6ab 2 u 3 



+ T, \ 6a(l + a) + (8 + 24a - 3a 2 + 8a 3 )u c 

(1 + a)f' 2 {u c ) L 

- 18(a(l + a) 2 )u 2 c + 15a 2 (l + a)ul }fc 4 + • • 



0. 
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